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In our previous study, we found two rbcL types (JN and JS types) in Japanese Conocephalum japonicum 
Grolle. In this study, we collected plant materials from 20 localities in China, Korea and Taiwan. We 
determined 1304 bps of each nucleotide rbcL sequence and found two rbcL types in Asian C. japonicum. 
One of them was a new rbcL type (CS type), which differs from the JN and JS types by 10 and 8 
nucleotides, respectively. The CS and JS types were found both in Mainland China and Taiwan. We could 
not find any morphological differences between the two types. In Chitou, Taiwan, the two rbcL types 
were growing together within a small research area measuring 100 m x 200 m. We performed an 
allozyme analysis on this population and observed that the CS type had Est-4a and Tpi-la alleles; where¬ 
as the JS type had Est-4b and Tpi-lb. It was strongly suggested that these two rbcL types might be repro- 
ductively isolated cryptic species of C. japonicum. 
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Bryophytes are typically the smallest land plants 
having a very simple morphology. However, most 
species of bryophytes have been recognized and 
described mostly by their limited morphological 
information. Morphological species of bryophytes 
often have a broad geographical distribution span¬ 
ning two or more continents (Schofield & Crum 
1972). Even though bryophytes can be easily dis¬ 
persed to remote sites through minute spores, such 
ranges seem too wide for the disrtribution of a sin¬ 
gle biological species. Therefore, we can expect 
that several cryptic species, which are genetically 


and ecologically well-differentiated from each other 
but morphologically hardly distinguished, might 
be contained within a single morphological species. 
In fact, cryptic speciation has been demonstrated in 
some broadly distributed bryophyte species such 
as the thallose liverworts ( Conocephalum conicum 
Dumort., Marchantia polymorpha L., Reboulia 
hemisphaeriea Raddi, Riccia dictyospom Howe, 
Aneura pinguis Dumort., Pellia epiphylla Corda), 
and a leafy liverwort ( Porella platyphylla Pfeiff) 
as well as several mosses ( Climacium americanum 
Brid., Leucobryum glaucum Angstr., Plagiomunium 
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cuspidatum T.J.Kop., Mielichhoferia elongata Nees 
& Hornsch. and Fontinalis antipyretica Hedw). 
(reviewed by Shaw 2001). These cryptic species 
have mainly been determined through allozyme 
analysis. 

Among these examples, the first known case of 
cryptic speciation in bryophtyes was demonstrated 
using the thallose liverwort Conocephalum 
conicum. Odrzykoski & Szweykowski (1991), 
Akiyama & Hiraoka (1994) and Kim et al. (1996) 
examined allozyme variation from throughout the 
geographical distribution range of the species, and 
found six discrete groups in total. The genetic dif¬ 
ferentiation among these six groups as measured by 
Nei's genetic distance was as large as those com¬ 
monly reported between different vascular plant 
species and much larger than those between con- 
specific populations of a higher plant species. They 
found several mixed populations of two different 
allozyme types of C. conicum , but even in these 
sympatrical populations, recombinant gametophytes 
were not detected (Szweykowski et al. 1981, 
Odrzykoski 1987). They concluded that these six 
allozyme groups were biologically independent 
cryptic species. 

In our previous study (Miwa et al. 2003), we 
sought to find cryptic species in Conocephalum 
japonicum Grolle using rbcL sequence variations 
rather than allozyme analysis as a preliminary indi¬ 
cator. In angiosperms, rbcL variations have been 
considered to be useful only for inter-generic or 
inter-familial level of phylogenetic analyses (Chase 
et al. 1993). In various ferns, however, rbcL varia¬ 
tion has been shown to be a useful tool for finding 
cryptic species (Murakami et al. 1998a, b, Yatabe et 
al. 1998, 1999, 2001, Kato et al. 2001, Masuyama et 
al. 2002). We considered it highly possible that 
rbcL variation would also be useful for finding 
cryptic species of bryophytes. 

Actually in our previous study (Miwa et al. 
2003), two distinct types of rbcL sequences (JN 
and JS), differing by 6 nucleotides, were found in 
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Japanese Conocephalum japonicum. Their geo¬ 
graphical distribution areas in Japan were clearly 
separated by the band of latitude 42-43 °N through 
Hokkaido. The geographical boundary of the two 
rbcL types is nearly identical to the northern limit of 
the distribution range of Fagus crenata in Japan 
(Horikawa 1972). We found a mixed population of 
JN and JS types in Abuta, Hokkaido, Japan. At this 
site, the two rbcL types were growing side by side in 
2-3 meters distance. We also conducted allozyme 
analysis and detected polymorphisms in EST and 
TPI enzymes, and tight associations between the 
allozyme and rbcL variants were observed; the JS 
type were fixed for Est-3a and Tpi-la alleles, where¬ 
as the JN type for Est-3b and Tpi-lb. Individuals 
with genotypes of other combinations were not 
found in this population. Thus, we tentatively con¬ 
cluded that these two rbcL types of C. japonicum 
represented two reproductively isolated cryptic 
species. 

Conocephalum japonicum is widely distrib¬ 
uted in East Asia: Nepal, Japan, Kamchatka (Russia) 
and Taiwan mark the west, east, north, and south 
boundaries of its range (Kitagawa 1982). In the 
present study, we collected samples of C. japonicum 
in China, Korea and Taiwan in order to discover pos¬ 
sible new rbcL types and also to elucidate the whole 
distribution range of the two rbcL types found in 
Japan. 

Materials and Methods 

Plant materials 

We collected Conocephalum japonicum at 20 
remote localities in China, Korea and Taiwan (Table 
1). We generally collected a single individual from 
each locality. For one population at Chitou, Taiwan, 
we mapped all the recognizable individuals growing 
in a 100 m x 200 m transect, and then collected a 
small amount of each of 21 plants for rbcL sequenc¬ 
ing and enzyme electrophoretic analysis. All vouch¬ 
er specimens are deposited in the Herbarium of the 
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Graduate School of Science, Kyoto University 

(kyo). 

DNA sequencing and molecular phylogenetic analy¬ 
sis 

Total DNAs were extracted using the Nucleon 
PhytePure plant and fungal DNA extraction kit 
(Amersham), or alternatively the Plant DNeasy Kit 
(Qiagen). Partial rbcL gene segments were amplified 
by PCR using Ready-To-Go PCR Beads 
(Amersham). We used two sets of primers: 1-1 
(ATGTCACCACAAACAGAGACTAAAGC), 2R 
(CTTCTGCTACAAATAAGAATCGATCTCTC- 
CA), N2-1 (TGAAAACGTGAATTCCCAAC- 
CGTTTATGCG), NN3-2 (GCAGCAGCTAGTTC- 
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CGGGCTCCA) for the PCR (Hasebe et al. 1994). 
A typical PCR amplification included an initial 
denature (5 min, 94°C) followed by 35 cycles with 
a 1 min denature at 94 °C, 1 min annealing at 50- 
55°C, 2 min 30 sec synthesis at 72 °C, and a final 
step of synthesis for 6 min at 72 °C. The PCR prod¬ 
ucts were purified using a QIA quick PCR 
Purification Kit (Qiagen). The amplified fragments 
were sequenced with a Big Dye terminator cycle 
sequencing kit (Applied Biosystems) using the 
above mentioned primers, and run on an Applied 
Biosystems Model 377 automated DNA sequencer 
(Applied Biosystems). The obtained sequences were 
handled by Sequence Navigator software (Applied 
Biosystems). A molecular phylogenetic tree was 
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Table 1. Voucher information of the plant materials from China, Korea and Taiwan. 
Abbreviation of collectors: XC, X. Cheng; TH, T.-W. Hsu; HM, H. Miwa; JS, J. 
Suhara. 


Locality 

Collectors 

Specimen No. 

C.japonicum (CS type) 

1. Kunming, Yunnan, China [2000 m] 

JS and XC 

HM4002 

2. Wuding, Yunnan, China [2000 m] 

JSandXC 

HM1272 

3. Daguan, Yunnan, China [1300 m] 

JSandXC 

HM3305 

4. Shilin, Yunnan, China [1500 m] 

JS and XC 

HM4014 

5. Malipo, Yunnan, China [750 m] 

JS and XC 

HM4061 

6. Xichou, Yunnan, China [1500 m] 

JSandXC 

HM4062 

7. Hekou, Yunnan, China [1150 m] 

JSandXC 

HM4006 

8. Taipingshan, Ilan, Taiwan [890 m] 

HM and TH 

HM3351 

9. Taipingshan, Ilan, Taiwan [1840 m] 

HM and TH 

HM3356 

C. japonicum (CS/'JS type) 

10. Chitou, Nantou, Taiwan [ca.l 150 m] 

HM and TH 

HM3209C1-21 

C. japonicum (JS type) 

11. Kochiu, Yunnan, China [1490 m] 

JSandXC 

HM4013 

12. Manhaozhen, Yunnan, China [540 m] 

JSandXC 

HM4010 

13. Adebo, Yunnan, China [1335 m] 

JSandXC 

HM4011 

14. Jinpin, Yunnan, China [1880 m] 

JSandXC 

HM4012 

15. Pinbian, Yunnan, China [1250 m] 

JSandXC 

HM4005 

16. Tianwangmiao, Taipei, Taiwan [470 m] 

HM and TH 

HM3361 

17. Mucha, Taipei, Taiwan [300 m] 

HM and TH 

HM3391 

18. FuanHuag, Nantou, Taiwan [865 m] 

HM and TH 

HM3381 

19. Cheju, Cheju-isl., Korea [600 m] 

HM 

HM3111 

20. Cheju, Cheju-isl., Korea [720 m] 

HM 

HM3113 
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constructed using the CLUSTALX program, version 
1.8 (Thompson et al. 1997). 

Allozyme analysis 

About 100 mg of fresh thallus tissue was homoge¬ 
nized in a 1 ml extraction buffer (0.2% 2-mercap- 
toethanol, 2% polyvinyl-pyrrolidone, 0.1 mM Tris- 
HC1, 1 mM EDTA4Na, 10 mM KC1, 10 mM 
MgC12, pH = 7.5). After centrifuging the 
homogenates at 10,000 rpm for 10 minutes at 4°C, 
20 pi of the supernatant was used for electrophore¬ 
sis for detecting each enzyme. Vertical polyacry¬ 
lamide slab gel electrophoresis using a discontinu¬ 
ous buffer system (Shiraishi 1988) was conducted 
for genotyping all plant samples collected from the 
Chitou population. Sixteen enzyme systems were 
tested for activity, and scorable bands were obtained 
for the following 10 enzymes: fluorescent esterase 
(EST E.C.: 3.1.1.1), triosephosphate isomerase (TPI 
E.C.:5.3.1.1), 6-phosphoglconate dehydrogenase 
(6PG E.C.: 1.1.1.44), shikimate dehydrogenase (SKD 
E.C.: 1.1.1.25), glutamate dehydrogenase (GDH 
E.C.:1.4.1.2), diaphorase (DIA E.C.:1.6.4.3), phos- 
phoglucoisomerase (PGI E.C.:5.3.1.9), glutamate 
oxaloacetate transaminase (GOT E.C.:2.6.1.1), acid 
phosphatase (ACP E.C.:3.1.3.2), and superoxid- 
edesmutase (SOD E.C.: 1.15.1.1). The loci were 
numbered in relation to their mobility: the longest 
distance zone was isozyme 1, and shorter distance 
zones were isozymes 2, 3 and so on. The allele 
labels alphabetically denote relative mobility, the a- 
alleles of each locus have the longest distance, fol¬ 
lowed by b, c and so on. 

Results 

We determined 1,304 nucleotide sequences of the 
rbcL gene for all the collected samples of Cono- 
cephalum japonicum. Two distinct rbcL sequences 
were found in Asian C. japonicum (Fig. 1). One of 
them was a new rbcL type, which differs from JN 
and JS types by 10 and 8 nucleotides, respectively. 


We named the new rbcL type as CS (China South) 
type. The CS type was found in the northern part of 
Yunnan-province, China and at higher elevations 
of Taiwan. The JS type was found in Korea, 
Yunnan-province, China and Taiwan. The geo¬ 
graphical distribution of CS and JS types is shown in 
Fig. 2, along with the distribution data of the JS and 
JN types in our previous work (Miwa et al. 2003). 
The JN type was not found among the samples in 
Mainland China and Taiwan. 

In Chitou, Taiwan, we found a mixed popula¬ 
tion of JS and CS types (16 individuals of JS type 
and 5 individuals of CS type). In this site, the two 
rbcL types were growing side by side within a 
small research area measuring 100 m x 200 m (Fig. 
3). We performed an allozyme analysis of the 21 
samples collected from the Chitou population. 
Polymorphisms were detected only in EST and 
TPI. In addition, the JS type always had Est-4b 
and Tpi-lb alleles, whereas the CS type always 
had Est-4a and Tpi-la, and thus tight associations 
between allozyme and rbcL variations were 
observed (Fig. 4). The appearances of the three 
rbcL types are shown in Fig. 5. Following the study 
on Conocephalum conicum by Odrzykoski & 
Szweykowski (1991) and Akiyama & Hiraoka 
(1994), we observed morphological features of cells 
and air chambers. We also compared length, width, 
thickness of thalli as well as their color among JN, 
JS and CS types. We could not find any significant 
morphological differences among the three rbcL 
types. 

We made a molecular phylogenetic analysis 
based on the obtained rbcL sequences of Cono¬ 
cephalum japonicum together with those of other 
bryophytes including C. conicum from the DNA 
database (Fig. 6). The rbcL sequences of C. japon¬ 
icum and C. conicum made a separate clade on the 
obtained molecular tree. Among the three rbcL 
types of C. japonicum, the JS and CS types were 
shown to be more closely related to each other than 
either is to the JN type. 
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JN) 

JS) 

CS) 


-GGATTCA AAGCTGGTGT TAAAGATTAT AGATTAAOIT ATTACACTCC 80 


JN) 

JS) 

CS) 


GGATTATGAA ACTAAAGAGA CAGATATTTT AGCAGCATTT AGAATGACTC CTCAGCCTGG GGTACCAGCA GAAGAAGCAG 160 


JN) 

JS) 

CS) 


GAGCAGCAGT TGCTGCTGAA TCTTCTACCG GTACATGGAC TACAGTTTGG ACTGATGGTC TTACTAATCT TGATCGTTAT 240 


JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

csl 

JN) 

JS) 

CS)_ 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 

JN) 

JS) 

CS) 


AAAGGTCGAT GCTATGGTAT TGACCCTGTT BCTGGAGAAG AAAATCAATA TATTGCTTAT GTAGCTTATC CTTTAGATTT 320 


ATTTGAAGAA GGATCTGTTA CAAATATGTT TACTTCTATT GTAGGTAATG TATTTGGATT TAAAGCTTTA AGAGCATTAC 400 


GTCTTGAAGA TTTACGAATT CCTCCAGCTT ATACAAAAAC TTTCCAAGGT CCTCCTCATG GTATTCAAGT TGAAAGAGAT 480 


AAATTAAACA AATATGGTCG TCCTTTATTA GGATGTACTA TTAAACCAAA ATTAGGTTTA TCTGCTAAAA ATTATGGTAG 560 


AGCTGTATAT GAATGTCTTC GTGGTGGACT TGACTTTACT AAAGATGATG AAAATGTAAA TTCTCAACCA TTTATGCGTT 640 


GGAGAGATCG TTTCTTATTT GTAGCAGAAG CTCT 



)CA GAAACTGGAG AAATCAAAGG ACATTATTTA 720 


AATGCTACTG CAGGTACATC TGAAGAAATG TTAAAAAGAG CAGCATGTGC TAGAGAGTTA GGTGTACCAA TTGTTATGCA 800 


TGACTATTTA ACTGGTGGTT TTACTOCAAA TACTAGTTTA TCTCATTATT GCCGTGATAA TGGTTTACTT CTTCATATTC 880 


ACCGTGCAAT GCATGCAGTT ATTGATAGAC AAAAAAATCA TGGTATCCAT TTCCGTGTAT TAGCTAAAGC TTTACGTTTA 960 


TCTGGTGGAG ACCATATTCA TGCTGGTACT GTTGTAGGTA AACTTGAAGG AGACCGT) 



iGTAGA 1040 


TTTACTTCGT GATGATTATA TCGAAAAAGA TAGAAGTCGT GGTATTTATT TTACACAAGA TTGGGTTTCT TTACCAGGTG 1120 


TTCTACCTGT AGCGTCTGGA GGAATTCATG TTTGGCATAT GCCTGCTTTA ACTGAAATTT TTGGAGATGA CrCTGTTTTA 1200 

C 


CAATTTGGTG GTGGAACT7T AGGC G CACCTGGTGC AGTTGCTAAC CGAGTTGCTT TAGAAGCGTG 1280 



TGTACAAGCA CGTAATGAAG GTCGTGATCT TGCTCGTCJNA GGTAATGAAA TTATTCGTGA 


1360 


Fig. 1 . Nucleotide sequences of the rbcL gene of three types in Asian Conocephalum japonicum. The variable sites are enclosed in 
boxes. The sequence data have been deposited in the DDBJ DNA database with the accession number JN type: AB046694, JS 
type: AB046695 and CS type: AB046719. 
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Fig. 2. The geographical distribution of the CS type (★), JS type (A) and JN type (♦) of Conocephalum japonicum in Asia. Details 
for Japanese populations are shown in Miwa et al. (2003). 



Fig. 3. Map showing the distribution of the CS (★) and JS type 
(A), of Conocephalum japonicum in the Chitou population, 
Taiwan. 


Discussion 

The 8-10 substitutions among the three rbcL types 
of Conocephalum japonicum may represent a large 
genetic distance considering the slow evolutionary 
rate of the rbcL gene of higher plants (Chase et al. 



C.S-lyl«: JSiypc CS typo JS ty|«' 


Fig. 4. The allozyme phenograms (EST and TPI) of the CS and 
JS types of Conocephalum japonicum in the Chitou popula¬ 
tion, Taiwan. 


1993, Yatabe et al. 1999). In order to test the 
hypothesis that these three rbcL types are repro- 
ductively isolated cryptic species, it is important 
to examine a sympatrical population where the two 
or more rbcL types were found growing together, as 
we did for the JN and JS types using their mixed 
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Fig. 5. General appearances of Conocephalum japonicum. A, JN type from Hokkaido, Japan. B, JS type from Kanagawa, Japan. C, 
CS type from Ilan, Taiwan. Bar, 1 cm. 


population in Abuta, Hokkaido, Japan (Miwa et al. 
2003). In this study, the allozyme analysis for the 
Chitou population showed that the JS and CS types, 
both recognized by chloroplast gene differentia¬ 
tion, were genetically fixed to different alleles in the 
Est-4 and Tpi-1 loci. Thus, genetic differentiation in 
two nuclear genes between the two rbcL types was 
clearly shown in this study as in the case of JN 
and JS types. Their results suggested that the three 
rbcL types (JN, JS and CS) are different biological 
species, because the results based on two indepen¬ 
dent gene markers, cpDNA {rbcL) and allozyme 
(nuclear), showed strong correlation. 


At this moment, we can not exclude the pos¬ 
sibility that all individuals of the Abuta and Chitou 
populations might be interfertile because Cono¬ 
cephalum japonicum produces numerous gemmae, 
and can propagate by an asexual process as well. We 
need to perform artificial crossing experiments in 
order to clearly show the existence of reproductive 
isolation among the three rbcL types. Anyway, we 
have found three genetically well-differentiated 
rbcL types in C. japonicum. 

Shaw (2001) noted that most of the cryptic 
species documented within widespread liverworts 
are genetically monomorphic or nearly so, in spite of 
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Fig. 6. Neighbor-joining tree based on the nucleotide sequences of the rbcL gene. Marchantia polymorpha (DNA database accession 
No X04465) and Dumortiera hirsuta (No U87068) were used as outgroups. We also added samples of Conocephalum conicum 
from five distant localities: WA, U.S.A. (No AB056154), Kanagawa, Japan (No AB046697), Hokkaido, Japan (No AB046696), 
NC, U.S.A. (No U87066), IL, U.S.A. (No U87067) for comparison. Numbers above the branches are the bootstrap values (%). 
Bootstrap analyses were accomplished using 1000 bootstrap replicates. 


high levels of genetic difference among them. He 
discussed how the origin of these cryptic species 
might be related to highly inbred or clonal repro¬ 
duction. The genetic situation found in Cono¬ 
cephalum japonicum is in fact the same as those 
generally found in other widespread liverworts. 
Capacity for clonal reproduction of C. japonicum 
might be related to its cryptic speciation. 

This study also suggested that rbcL sequence 
variation is useful in searching for cryptic species of 
bryophytes. Many cryptic species of bryophytes 
have been found using allozyme analyses (see Shaw 
2001). However, cryptic species-specific alleles 
were rarely found, and it is usual that allozyme 
variations are observed in each cryptic species. To 


date, cryptic species recognized by allozyme analy¬ 
ses have generally been identified as clusters in 
UPGMA dendrograms based on Nei's genetic dis¬ 
tances (Dewey 1989, Odrzykoski & Szweykowski 
1991, Akiyama & Hiraoka 1994). This indicates 
that bryophyte cryptic species were first recognized 
by differentiation in allele frequencies of several 
independent allozyme loci. Cryptic species-specif¬ 
ic allozyme bands can always be found and identi¬ 
fied after clusters of their population have been 
recognized. In other words, allozyme analyses usu¬ 
ally help in identifying only populations, even when 
the analysis is a search for cryptic species. 
Therefore, if plants of different cryptic species are 
growing sympatrically, as in the case of JS and CS 
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types in the Chitou population, it may be difficult to 
find cryptic species using allozyme analyses. 

In such cases, we expect rbcL sequence analy¬ 
sis to be useful because it is always applicable on the 
individual level. The rbcL genes are encoded by 
chloroplast DNA, which is usually treated as mater¬ 
nally inheriting without recombination. Maternal 
inheritance of the plastid genome was also docu¬ 
mented in one species of bryophytes ( Anthoceros 
punctatus L.) by Izumi & Ono (1999). In this case, 
effective population size for maternally inherited 
genes is expected to be one quarter of that for the 
biparentally inherited nuclear genes such as 
allozyme loci (Moore 1995). Therefore, genes 
encoded by cpDNA are more easily fixed in each 
Mendel population by random genetic drift. This 
means that their variations, if available, tend to 
coincide with the boundaries of other biological 
species even in mixed populations. 

Moreover, samplings are much easier in rbcL 
sequencing than in allozyme analysis, because even 
silica gel dried samples are sufficient for the former 
method, whereas fresh materials are really necessary 
for the latter. For allozyme analyses, standard plant 
samples are also necessary for identifying allozyme 
bands in different electrophoresis, but it is difficult 
and laborious to keep bryophyte samples alive for 
long periods. Such procedures are completely 
unnecessary for rbcL analyses. 

These advantages of cpDNA sequence data 
suggest rbcL sequences to be more useful than 
allozyme data as a preliminary indicator of repro- 
ductively isolated assemblages or biological species 
(Yatabe et al. 2002). In fact, rbcL variations have 
been shown to be a good key for finding cryptic 
species in various fern groups. The present study 
suggests that rbcL variations are also useful for 
identifying cryptic species of bryophytes. We con¬ 
sider it necessary to examine actual circumstances of 
rbcL variation among the cryptic species recog¬ 
nized by allozyme analyses in previous cases. In 
addition, usefulness of rbcL variation in searching 
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for cryptic species will be also elucidaced if we 
find a new one in Conocephalum conicum, in which 
cryptic species have been well sought by allozyme 
analyses. It is indeed, our next step to establish the 
new, more efficient method for recognizing bio¬ 
logical species in bryophtytes. 
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